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Abstract 

In this communication, the band-gap structures of n-CF piezoelectric superlat- 
tices have been calculated using the transfer-matrix-method, the self-similarity 
behavior and recovery rule have been systematically analyzed. Consistent with 
the rigorous mathematical proof by Hu et al.[A. Hu et al. Phys. Rev. B. 48, 
829 (1993)], we find that the n-CF sequences with 2 < n < 4 are identified 
as quasiperiodic. The imaginary wave numbers are characterized by the self- 
similar spectrum, their major peaks can all be properly indexed. In addition, we 
find that the n = 5 sequence belongs to a critical case which lies at the border 
between quasiperiodic to non-quasiperiodic structures. The frequency range of 
self-similarity pattern approaches to zero and a unique indexing of imaginary 
wave numbers becomes impossible. Our study offers the information on the 
critical 5-CF superlattice which was not available before. The classification of 
band-gap structures and the scaling laws around fixed points are also given. 
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1. INTRODUCTION 

From a symmetrical point of view, the solids in nature are classified as 
crystals, quasicrystals, and noncrystalsflj. The crystal structure is character- 
ized by a periodic arrangement of atoms in space, the waves propagating in 
crystals have a Bloch wave form. According to the Bloch theorem all states 
are extended in space and eigenenergies form continuous band structures. For 
fully disordered noncrystal, it is believed that all states are localized in space 
and energy spectra are singular continuous. Quasicrystal is a unique type of 
structure which lies at the boundary between translation invariant crystals and 
random glassy materials Q. Though lacking long-range translational symmetry, 
it possesses a certain orientational order. To understand the interim situation 
between crystals and noncrystals, much research works have been carried out on 
quasicrystals, especially since the experimental discovery of quasicrystal phase 
in an Al-Mn alloy with icosahedral symmetry The simplest structural model 
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describing quasicrystals is the Fibonacci lattice model. It has received a great 
deal of attention^ S, S 0, i, II, E3, EI El El, E! El since it contains the basic 



ingredients of quasicrystals and is relatively easy to deal with. The Fibonacci 
lattice model is famous for its Cantor-set spectrum, the self-similar spectrum is 
identified as critical eigenstates which is neither extended as in periodic system, 
nor localized as in disordered onesQ- 

To explore the spectral evolution from the translational invariant crystals to 
fully disordered materials, a structural model which is tunable between the two 
extreme cases are highly desirable. A simple Fibonacci lattice model is certainly 
a good candidate since its symmetry lies in between, but it is not good enough 
since its order is not adjustable. In this respect, the generalized Fibonacci 
lattice model with n-component (n-CF) offers an excellent choice and its order 
can be tuned by the parameter n[lj|. The n-CF lattice can be obtained by 
repeated application of the concurrent substitution rules TA\ = A\A n , TAi = 
Ai-i(i — 2, . . . , n), where Ax, A2, ■ . . , A n are the n-components and T is the 
substitution operator. The Ith generation sequence of n-CF lattice is defined 
as Si — T l A\ and the generalized n-CF lattice corresponds to Sao. The first 
few generation sequences are given by So = A\, S% = AiA n , S2 = AiA n A n ^x, 
and S„ = AiAnAn-i . . . A3A2A1. In general, Si = Si-i + Si- n (l > n). 
The generalized n-CF lattice recovers the periodic lattice if n is set to n = 1, it 
also naturally reproduces the standard Fibonacci lattice when n = 2. What's 
unique about the generalized n-CF lattice is that it can be tuned smoothly 
from periodic, to quasiperiodic, and then finally to non-quasiperiodic lattices as 
n changes. As shown by Hu et ai.[lf3], the information on the quasiperiodicity 
of a n-CF lattice is encoded in the substitution matrix 



T = 



( 1 ... 1 \ 
10 ... 
1 ... 



V ... 1 oJ n 



(1) 



whether the n-CF lattice is quasiperiodic or not is closely associated with the 
Pisot number of characteristic polynomial equation A™ — A™ -1 — 1 = 0. A Pisot 
number is present if the characteristic polynomial equation has only one root 
I < A ™' < 2 and all its conjugate roots have their norms less than 1. Following 
the Bombier- Taylor theorem, the existence of Pisot number guarantees that the 
n-CF lattice is quasiperiodic and can be generated by the cut and projection 



method 17]. Hu et o/.[16j showed that quasiperiodic lattice exists only if 2 < 
n < 5; n = 1 corresponds to a periodic lattice while lattices with n > 5, though 
still ordered, is no longer quasiperiodic. In fact, a close inspection of the root 
structure suggests that n = 5 lattice should be classified as a critical case since 
there is a pair of conjugate roots having their norms equal 1, which neither 
belongs to the quasiperiodic case where all conjugate roots have their norms 
less than 1, nor belongs to a non-quasiperiodic case where some of the conjugate 
roots have their norms larger than 1. The Pisot numbers are A = 1.61803, 



2 



A* 3) = 1.46557, A 



(4) _ 



1.38028, and A, 



(5) _ 



sQl, \m [io, m m, M, H3, 



, - 1.32472 for n = 2, 3, 4, 5, respectively. 
Above discussion suggests that n-CF lattice covers a wide range of lat- 
tice types from periodic, quasiperiodic, critical, and non-quasiperiodic lat- 
tices, it offers a natural platform to study the spectra evolution with struc- 
tural ordering. Though extensive investigations have been carried out on the 
electronic di,i, 0,11, vibrational d [13, [ul, and dielectric [H, [H, E3, El prop- 
erties of quasiperiodic structures as a vehicle to study the evolution process from 
ordered periodic structures to disordered solids, the most previous studies con- 
centrated on the standard Fibonacci lattices [H, [jl, El [7 , |j , [lO . 
and the generalized 2-component Fibonacci lattices 
the generalized n-CF lattices are mostly discussed in the context of structural 
propertics[25]. Furthermore, the previous studies mostly dealt with the single 
degree problems, and studies on the mode-coupling problems with energy trans- 
fer between different degrees of freedom are only a few 2f| 27, 28|. Thus, we 
have carried out a comprehensive study in this communication on the band-gap 
structures of n-CF piezoelectric superlattices. The transfer-matrix-method is 
explored so that the band-gap structures can be most easily visualized through 
the imaginary wave number. Our study show that the band-gap structures can 
be classified into a series of big cluster in the frequency spectrum. For each n- 
CF superlattice, the self-similar spectrum and recovery rule are present within 
each generation and across different generations I. Consistent with the rigor- 
ous mathematical proof, we find that the n-CF superlattices with 2 < n < 4 
can be identified as quasiperiodic. The band-gaps are characterized by the 
self-similarity behavior and all the gaps can be properly indexed; The n = 5 
superlattice belongs to a critical case which lies at the border between quasi- 
periodic to non-quasiperiodic structures. The frequency range of self-similar 
spectrum approaches zero, a unique indexing of band-gaps becomes impossible. 



2. MODEL 

In this communication, we use superlattices made of piezoelectric LiNb03 
compound as our model systems [2a], To distinguish different components of 
n-CF superlattices [l6j, each component is made of a pair of positively and neg- 
atively polarized ferroelectric domains of varying lengths (L\, L' l _). For sim- 
plicity and also for setting up length scale, the lengths of positively polarized 
domains of all components are taken to be the same L+, while the lengths of 
negatively polarized domains are used to differentiate the different components. 
A typical superlattice for a given n-CF sequence and for a given generation I 
is illustrated in Fig. 1. The experimental setting and material parameters are 



described in detail in the references. 27 



To calculate the band-gap structures of n-CF lattices, we use the sequence of 
Ith generation as a supercell to form a periodic lattice. In this case, the lattice 
constant ai = J27=i wnere F\ is the number of ith component in Ith gen- 
eration sequence and L l = ttL 1 / L + is the reduced length of ith component [271]. 
The corresponding reciprocal lattice vectors are K mi = (27r/a;)m;s with mi 
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denoting an integer, the reduced wave number —n/ai < k < +n/ai. Thus, 
the essential information on the Ith generation of n-CF sequence is included in 
its Fourier components 6 (mi). The band-gap opens at the reduced frequency 
uj « K mi if 8 (m{) / 0. For the special case where L\ = it and L l _/iv can all 
be written as rational numbers, L % — tt[1 + V_ja % J\ with the greatest common 
divisor gcd(a l _ , b l _) = 1, the reciprocal lattice vectors can also be written as 
rational numbers 

= 2m; x lcm(aL,ag.,---,a") 

mi ~ lcm(aL,a^,---,a") + &-/ a -) [ ' 

with lcm(a?_, • • ■ , a™ ) denoting the least common multiple. It is easy to 
shown that 9 (mi) vanishes when mi = integer x lcm(aL, a_, ■ • • , a") -J^™ + 
b t _/a i _). In this case, the band-gap structures are divided into big clusters, the 
self-similarity and scaling behavior in each cluster can be used to check whether 
the given superlattice is quasiperiodic or not. 



3. NUMERICAL RESULTS AND DISCUSSIONS 

To study the spectral evolution with generation I for n-CF superlattices, usu- 
ally the furcation pattern of eigenmodes is traced to analyze the self-similarity 
property and to search for the fixed point. When system involves continuum 
eigenmode spectra, it is better to track the furcation pattern of band-gap in- 
stead of the eignemodes since the discrete band-gap structures encode the essen- 
tial information of n-CF sequences. As the band-gap represents the forbidden 
frequencies of eigenmodes, the wave propagating in this frequency range is at- 
tenuated at an imaginary wave number Imfc. For polaritons which involve the 
coupling between electromagnetic wave and acoustic wave, the effective dielec- 
tric function in the vicinity of a band-gap always diverges at low band-gap edge 
while continuously approaches to zero at high band-gap edge. This suggests 
that the asymmetrical peaks symbolizing imaginary wave numbers can be used 
to visualize the band-gap structures. In the following, we shall systematically 
analyze the band-gap spectra of n-CF superlattices as functions of generation I 
and component n. 

Since the 2-CF sequence corresponds to the standard Fibonacci lattices 
which have been extensively studied before[28j], we concentrate on 3 < n < 5 
cases below. In choosing the components of n-CF sequence, we follow the rule 
so that L 1 > L n > • • • > L 2 16]. In Fig. 2, the imaginary wave numbers 

Imfc as functions of frequency u) are presented for the 9th, 12th, and 16th gen- 
erations of 3-CF superlattices. The three components are chosen as L x _ = tt, 
L 2 _ = l7r/3, and L z _ = 2ir/3. The overall patterns of the band-gap spectra 
look quite similar among different generations and the positions of major band- 
gaps are essentially independent of generations. Within the same big cluster, 
the spectra demonstrate a certain symmetry with respect to the center of big 
cluster(o) = 3). In addition, the big cluster is further subdivided into three 
smaller clusters and each subcluster has its spectral center. The spectral detail 
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increases with generation I for a given n, thus the self-similarity pattern and 
scaling law can be analyzed when / is large enough. 

To check the quasiperiodicity of 3-CF superlattice, we follow Hu et al.'s 
procedure[3| and identify the major peaks of n-CF superlattice. According to 
the cut and projection method (l7j. the reciprocal lattice vectors are given by 



where D ~ Yli=i L l i]i is an average lattice constant and hiS are any integers. 
r/i is the relative weight of ith component and satisfies l/rj n = rj n /rj n ^i = 
■ ■ ■ = 113/112 = Aq , and Aq is the Pisot number of substitution matrix of n-CF 
sequence. For quasiperiodic superlattices of n-CF type (2 < n < 5), D converges 
quickly for large generation I. D does not converge for aperiodic superlattices 
when n > 5. Since the band-gap takes place approximately at w ~ K mi , one 
can properly identify the various major peaks by using K mi = k[hi, h,2, ■ • • , h n ]. 
In Fig. 2c, the major peaks of 3-CF superlattice of large generation are labeled 
in this way. To investigate the self-similarity property of 3-CF superlattices, 
the spectra of imaginary wave numbers are presented in Fig. 3 for the 9th and 
16th generations. With regarding to both the positions and magnitudes of the 
peaks, there is a clear one-to-one correspondence between them in the vicinity 
of the cluster center. The frequency range of the self-similar patterns is quite 
broad extending from u> — 2.29 to Q = 3.71. The dynamical property of 3-CF 
superlattices is recovered after 7-generation and the scaling parameter is (A ) 7 . 

Since the peak labeling of band-gaps is one of the essential feature of 
quasiperiodic superlattices, it is of interest to see how the situation goes for 
the 4-CF and 5-CF superlattices. In Fig. 4a, the spectrum for 12th genera- 
tion of 4-CF superlattice is presented, the four components are U_ — 47r/3, 
L 2 _ = br/3, L 3 _ = 2tt/3, L 4 _ = 3?r/3. While Fig. 4b is the similar spectrum for 
the 12th generation of 5-CF superlattice with component setting L\_ — 6n/3, 
L 2 _ = I71-/3, L 3 _ = 2tt/3, Li = 4tt/3, L 5 _ = 5tt/3. Labeling the band-gap peaks 
is a very tedious job for multi-component Fibonacci superlattices since more 
component means denser distribution of reciprocal lattice vectors. For the 4- 
CF superlattice, we are able to identify the major band-gap peaks with easy. 
As shown in Fig. 4a, all major peaks are uniquely labeled as it should be for 
the quasiperiodic lattice. The situation becomes more complex for 5-CF super- 
lattice since it lies at the border between quasiperiodic and non-quasiperiodic 
structures. We find that there exists three or more possible labelings for some 
major peaks, and their difference in wave number is less than 10~ 15 , smaller than 
the effective digits that double precision number can offer. This non-uniqueness 
seems peculiar to 5-CF superlattices, we believe that it may be related to its 
critical nature of the structure. 

Unlike the 3-CF superlattices whose band-gap spectra recovers every 7th 
generation around the fixed point at the cluster center, we find that the spectra 
for 4-CF superlattices recovers every 15th generation. Fig. 5 compares the imag- 
inary wave number spectra for both the 12th generation and expanded view of 




(3) 
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27th generation scaled by (A^) 15 . Due to one more component in composing 
the superlattice, the frequency range of self-similar patterns is significantly re- 
duced in comparison with that of 3-CF superlattices. But a clear one-to-one 
correspondence is still evident in the narrow range 2.865 < Oj < 3.135 both for 
peak positions and magnitudes. If we summarize the recovery rule get so far for 
the fixed point near the cluster center, following deductive formula seems rea- 
sonable Al = 2™ — 1. It predicts that the spectra recover every Al = 1, 3, 7, 15 for 
the periodic, standard Fibonacci [28j], 3-CF, and 4-CF superlattices. Of course, 
the above statement holds only for the fixed point near the spectral cluster 
center. As a matter of fact, there are many fixed points for a given dynamical 
system 29|, 3(|, and they do not share the same scaling parameter at all. 

If the deductive formula is right, one would expect to find the band-gap 
spectra of 5-CF superlattices recover every 31st generation, but we did not 
succeed. Though for different settings of components we are able to find an 
approximate fixed point around Q — 3 for the 3rd, 7th, and 21st generation 
recovery rule, none of them is perfect in the sense that either the fixed point 
is not strictly followed or the frequency range of self-similar spectra around the 
fixed point is too narrow. Thus, we are led to the conclusion that the critical 
5-CF superlattices are characterized by two important features: (1) the major 
peaks cannot be uniquely labeled; (2) the frequency range of self-similarity 
pattern shrinks to zero. 



4. CONCLUSION 

In this communication, the peak labeling and self-similarity patterns of imag- 
inary wave number spectra are used as essential tools to cross check the nature 
of n-CF superlattices. For n = 3 and n = 4-CF superlattices, major peaks can 
always be labeled uniquely. There exists finite frequency range of self-similarity 
pattern around the fixed point and the frequency range decreases as n increases; 
For the critical case of 5-CF superlattices, the unique labeling of major peaks be- 
comes impossible and the frequency range of self-similarity pattern approaches 
to zero. Our study substantiates the rigorous mathematical proof by Hu et 
a/.[T3], in particular, valuable information is supplemented on the critical 5-CF 
structures. 
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Figure 1: The components and piezoelectric superlattice of n-CF sequence, (a) All components 
have an identically positively polarized domain(L+ = n) plus a negatively polarized one with 
varying length, (b) An example of n-CF superlattice. 



Figure 2: The imaginary wave numbers as functions of frequency for 3-CF supcrlattices. The 
three components are L x _ = 7r, L 2 _ = 1tt/3, and = 2w/3. (a) For 9th generation; (b) 
For 12th generation; and (c) For 16th generation. For superlattice of large generation, major 
peaks can be labeled according to the projection and cut method. 



Figure 3: The self-similarity behavior around the fixed point of the cluster center in 3-CF 
supcrlattices. (a) For 9th generation; (b) For 16th generation, enlarged around u) = 3 and 

scaled by a factor (Aq 3 ') 7 . The other parameter settings are the same as Fig. 2. 



Figure 4: The indexing of major peaks of imaginary wave numbers, (a) For 4-CF. The four 
components are: V_ = 4vr/3, 1?_ = Itt/3, It = 2tt/3, L± = 3vr/3; (b) For 5-CF. The five 
components arc: V_ = 6?r/3, L 2 _ = Itt/3, L't = 2tt/3, L 4 _ = 4tt/3, L 6 _ = 5tt/3. For the 
critical case n = 5, the peak labeling is not unique anymore. 



Figure 5: The self-similarity behavior around the fixed point of the cluster center in 4-CF 
supcrlattices. (a) For 12th generation; (b) For 27th generation, enlarged around Co = 3 and 
scaled by a factor (Ay 4 ') 15 . The other parameter settings are the same as Fig. 4a. 
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This figure "Figure- l.png" is available in "png" format from: 



http://arxiv.org/ps/1012.0171v3 
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